Off shell behaviour of the in medium nucleon-nucleon cross 

section 

C. Fuchs^, Amand Faessler^, M. El-Shabshiry^ 

^ Institut fiir Theoretische Physik, Universitdt Tubingen, Tubingen, Germany 

'^Physics Department, Faculty of Science, Ain Shams University, Cairo, Egypt 
The properties of nucleon-nucleon scattering inside dense nuclear mat- 
ter are investigated. We use the relativistic Brueckner-Hartree-Fock model 
to determine on-shell and half off-shell in-medium transition amplitudes and 
cross sections. At finite densities the on-shell cross sections are generally sup- 
pressed. This reduction is, however, less pronounced than found in previous 
works. In the case that the outgoing momenta are allowed to be off energy 
shell the amplitudes show a strong variation with momentum. This descrip- 
tion allows to determine in-medium cross sections beyond the quasi-particle 
approximation accounting thereby for the finite width which nucleons acquire 
in the dense nuclear medium. For reasonable choices of the in-medium nu- 
clear spectral width, i.e. F < 40 MeV, the resulting total cross sections are, 
however, reduced by not more than about 25% compared to the on-shell val- 
ues. Off-shell effect are generally more pronounced at large nuclear matter 
densities. 
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I. INTRODUCTION 

One major topic of modern nuclear physics is the investigation of hadron properties inside 
a dense hadronic environment which exists, e.g. in the interior of neutron stars or is tran- 
siently created in energetic heavy ion collisions. In the latter case, the theoretical framework 
to describe the time evolution of heavy ion reactions is provided by kinetic transport theory. 
Starting from the quantum theory of strongly interacting Fermi systems, formulated within 
the framework of the Dyson-Schwinger hierarchy of non-equilibrium many-body Green func- 
tions one can derive semi-classical transport equations of a Boltzmann-Uehling-Uhlenbeck 
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(BUU) type [1,0. These transport equations describe successfully the time evolution of a 
non-equilibrated strongly interaction hadron gas. To mention only the essential steps of such 
a derivation there are: the truncation of the many-body hierarchy at the two-body level, 
a Wigner transformation of the density matrices with subsequent gradient expansion up to 
first order in % and the usage of the quasi-particle approximation (QPA) which neglects the 
finite decay width of the particles. The resulting BUU equation consists of two parts, a drift 
term which propagates the particles dressed by the surrounding medium in a self-consistent 
mean field, and the collisions term responsible for binary nucleon-nucleon scattering. In a 
consistent treatment both ingredients, namely the mean field and the binary cross sections 
should be treated on the same footing which means to base both of them on the same effec- 
tive interaction. Unfortunately, in most applications to heavy ion collisions this is not done. 
The self consistent mean field accounts for medium effects by its density dependence. For 
the cross section, on the other hand side, the free (vacuum) expressions are widely used in 
transport calculations. It has, however, been noticed that in particular at incident energies 
below the particle production thresholds, medium modifications of the cross sections can 
play an important role for the reaction dynamics in heavy ion collisions 

The relativistic (Dirac-) Brueckner approach |?|^T3[ provides a powerful tool to achieve 
such a consistent description. Starting from free nucleon-nucleon interactions, given in its 
modern form by one-boson-exchange potentials [1^ one treats the two-body correlations in 
dense nuclear matter in the ladder approximation of the Bethe-Salpeter equation. As a result 
the nuclear matter saturation properties are quite well described. This is achieved without 
the adjustment of additional parameters, as e.g. done in relativistic mean field models 
|T3|. On the level of the T- matrix approximation both ingredients for the BUU equation 
follow from the on-shell in- medium T-matrix (or G-matrix). The mean field is determined 
by the real part ReT of the T-matrix whereas the cross section o" ^ |Tp is connected to 
the imaginary part ImT via an optical theorem 0. Medium modifications arise due to 
the dressing of the quasi-particles and the existence of the Pauli operator which prevents 
the scattering of intermediate states in the Bethe-Salpeter equation (not final states) into 
occupied phase space areas. Both aspects are most pronounced at high densities and/or low 
momenta and lead to a suppression of the in-medium cross section compared to the free one. 
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There have aheady several studies been devoted to the in-medium NN scattering problem. 
The Tiibingen group ||^ and later on the Rostock group [|TB],|T^ derived in-medium cross 
sections within the non-relativistic Brueckner approach, in the latter case also at finite 
temperature. Relativistic calculations were performed in |[T8|-pO|. The most complete study 
of in-medium NN scattering within the Dirac-Brueckner approach was probably done by 
Li and Machleidt [^] who used the BONN potentials as bare interaction. Unfortunately 
the different approaches have led to partially contradictory results, in particular between 
relativistic |T^J2T|] and the non-relativistic calculations 0]. Therefore, in the first part of the 
paper we revisit the problem of on-shell scattering. The results are obtained with the BONN 
A potential. In large parts we find a good agreement with the previous investigations of Li 



and Machleidt [^. However, their treatment seems to overestimate the suppression of 
the in-medium cross section at low energies compared to the vacuum case. 

The second part of the paper is devoted to an additional aspect: Kinetic transport 
equations are essentially based on the quasi-particle approximation (QPA) which puts the 
nucleons on the mass-shell. The same holds in the medium for the dressed quasi-particles. If 
the imaginary part of the nucleonic self-energy is negligible (JmS <^ i?eE) the quasi-particle 
approximation as the zero-width limit (F oc JmS) for the nucleon spectral function appears 
to be justified 0. However, it is well known that the spectral widths of hadrons change in 
the medium. It has further been pointed out that in the medium also "stable" particles can 
obtain a non-zero width, depending on their collision rates (collisional broadening), see e.g. 
2^ . Following the work of Botermans and Malfiiet there have been several attempts 



to derive transport equations for non-equilibrated Fermi systems beyond the quasi-particle 
approximation, see and references therein. However, due to the complications which 

arise giving up the QPA these transport equations were never used in practical applications 
to heavy ion reactions but stayed more or less on the level of academical considerations. 
Just very recently there have been successful attempts to formulate generalised transport 
equations which can be handled in applications using testparticle methods [^,0. As the 
basic feature of these approaches energy and momentum of the testparticles are no more 
related by the mass-shell condition but according to their spectral distributions. Thus 
particles are propagated and also scattered off mass-shell. Hence, the knowledge of off- 
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shell scattering amplitudes becomes necessary. However, the behaviour of such amplitudes 
is presently unknown to large extent. The off-shell structure of the scattering amplitude 
determines in this context also the magnitude of non-local corrections to the Boltzmann 
equation which can be translated into non-local time and momentum shifts in the binary 
scattering process 



Here again the relativistic Brueckner approach provides the natural tool to determine 
in-medium off-shell scattering amplitudes in the ladder approximation. In the present work 
we extend the on-shell NN scattering to the half-off-shell case where the incoming parti- 
cles are still on their mass shells but the final states are generally off-shell. These matrix 
elements provide valuable information for future transport investigations beyond the QPA. 
We investigate the off-shell structure of the in-medium T-matrix, respectively the transi- 
tion amplitudes |Tp, and the resulting cross sections over a wide range of nuclear matter 
densities. 

The paper is now organised as follows: First (Sec. 11) we briefly sketch the basic features 
of the relativistic Brueckner approach. In Sec. HI the on-shell scattering problem in the 
medium is discussed. Neutron-neutron and proton-proton channels are considered separately 
and total and differential cross section are given. We also compare with results of other 



groups, mainly those of Li and Machleidt |2T|]. In Sec. IV we turn to the half off-shell 



case and discuss the structure of transition amplitudes and cross sections beyond the quasi- 
particle approximation and summarise in Sec. V. 

II. RELATIVISTIC BRUECKNER APPROACH 

In the relativistic Brueckner approach the Bethe-Salpeter equation is reduced to a three 
dimensional integral equation of the Lippmann-Schwinger type, the so called Thompson 
equation |2^. The Thompson propagator projects thereby the intermediate nucleons onto 



positive energy states and restricts the exchanged energy transfer by 6{k^) to zero. The 
Thompson equation is most easily solved in the two-nucleon c.m.-frame 

r rl^k A/f*2 (Kk T*) 

Tip, q. P)l.. = V(p. q) + / ^V(p. kj^jJ^^i^Tdc. q. P) (1) 
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where q = (qi — q2)/2 = (q^ — ^)/'^ is the relative momentum of the initial states and 
similar p, k are the relative momenta of the final and intermediate states, respectively. P = 
(qi + q2) is the cm. momentum. The starting energy in Eq. (|ID is fixed to = 2£'*(q). 
Sandwiching the one-boson-exchange potential V (|lD between in-medium spinors @) the 
matrix elements acquire a density dependence which is absent in non-relativistic treatments 
and which is believed to be the major reason for the much improved description of the 
nuclear saturation properties in the relativistic theory. The Pauli operator Q prevents 
intermediate states from scattering into forbidden phase space areas. 

Inside the medium the particles are dressed which leads to effective masses and the 
kinetic momentum 

M*(A;) = M + ReS,(A;) , k*" = + ReW{k). (2) 

Re and Im denote real and imaginary part since (above the Fermi surface) the self-energy 
is generally complex. Here we adopt the quasiparticle approximation, i.e., the Im [S] will 
be neglected in Eq. (|l]). This means that the decay width of the dressed nucleon state k 
to another state k' is set equal to zero, resulting in an infinite lifetime of this 'quasiparticle' 
state. Furthermore, the explicit momentum dependence of the self-energy which enters via 
a term kS^, proportional to the spatial component of the vector self-energy is small and 
can be dealt by introducing the reduced kinetic momentum k*'^ = k*^/{l + S^,) and the 
reduced effective mass M* = M*/(l + S^) 0]. Thus, the nucleons are given by plane waves 
which fulfil a quasi-free Dirac equation 

'^,k*^-M*]ux{k) = . (3) 

Using the normalisation of Ref. |^ the self-consistent positive-energy spinors of helicity A 
are defined as 

/ 



ux{k) 



E*{k) + M* 
\ 2M* 



1 

2A|k| 



XX (4) 



with xx being a Pauli spinor. The Dirac spinors depend on the effective mass and thus on 
the nuclear density. In the Thompson equation (|I]) and in the discussion below we deal with 
the rescaled quantities M* and k* but will omit this in the notation further on. 



To summarise the kinematics of the Thompson equation: 

1. The initial states are on-shell, i.e. = {-E'*(q), ±q} with E*{q) = \J M*'^ + = ^V^s*- 
The final states fulfil energy- momentum conservation = {^\fs*, ±p} and are off-shell as 
soon as |p| 7^ |q|. 

2. The determination of the off-shell matrix elements is perturbative in the sense that the 
quasiparticle approximation is applied to the Thompson equation, although T is generally 
complex for incident momenta above the Fermi surface which leads to a non-vanishing 
imaginary part of the self-energy JmS and correspondingly, an imaginary optical potential 

in. 



To determine the scalar and vector components Sq and S^, of the self-energy is a 
subtle problem. Here on-shell ambiguities arise form the projection onto positive energy 
states when the T-matrix is decomposed into Lorentz invariant amplitudes. This problem 
has been known for a long time and is still not completely resolved. In [0 we discussed 
the failure of previously used recipes |^,|TT| which lead to spurious contributions in the self- 
energy from the coupling to negative energy states, in particular spurious contributions from 
a pseudo-scalar one-pion-exchange which are not completely replaced by a pseudo-vector 
coupling. In Ref. this problem was extensively discussed and a method to minimise the 



on-shell ambiguities was proposed. Here we used the scheme of [|13[ where the Born term V 
and the remaining ladder kernel of the Thompson equation are treated separately. Thus we 
account properly for the pseudo-vector structure of the Born contributions \4 from vr and 
r] exchange contributions when the projection of the full T-matrix onto covariant amplitudes 
is performed. The remaining ladder kernel thereby is treated as pseudo-scalar. 

To solve the Thompson equation (|I|) in the c.m.-system we apply standard techniques 



which are in detail described by Erkelenz [^. After a partial wave projection onto the 
\JMLS >-states the integral reduces to a one-dimensional integral over the relative mo- 
mentum |k| and Eq. (|I]) decouples into three subsystems of integral equations for the 
uncoupled spin singlet, the uncoupled spin triplet and the coupled triplet states. Due to 
the antisymmetry of the two-fermion states we can restore the total isospin (1=0,1) of the 
two-nucleon system with the help of the selection rule (— = — 1 which means that 
matrix elements are already antisymmetrized. 
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The Pauli operator Q is replaced by an angle averaged Pauli operator Q. For non- 
vanishing cm. momenta the Fermi sphere is in the two-nucleon c.m.-frame deformed to a 
Fermi ellipsoid for which Q has to be evaluated [0,^ . We are solving the integral equations 
by the matrix inversion techniques of Haftel and Tabakin [^]. Real and imaginary parts 
of the T-matrix are calculated separately by the principal-value treatment given by Trefz 
et al. From there it is seen that ImT is essentially proportional to the angle averaged 
Pauli operator and thus it is strongly reduced for momenta below the Fermi surface due to 
Pauli blocking. Then positive-energy helicity T-matrix elements are constructed from the 
I J MLS >-scheme as described in |]T0| , |30| . From general symmetries it follows that for each 
total angular momentum J only six of the sixteen helicity matrix elements are independent 

ri=<+ + |r-^(p,g)| + +> , T2=<+ + |T-^(p,g)|--> 
T^=< + ~\T\p,q)\ + -> , T4=< + -|r-^(p,g)|- + > 
T^=<+ + \T\p,q)\ + -> ^T,=< + -\T\p,q)\ + +> (5) 

which in the on-shell case |p| = |q| further reduce to five independent matrix elements since 
then T5 = Tg. From the six independent amplitudes in the \ JMLS >-representation the six 
independent partial wave amplitudes (^) in the helicity representation (for 1=0,1 and real 
and imaginary part separately) are obtained as described in Ref. @]. Summation over J 
yields the full helicity matrix elements 

rO 7 -I- 1 "1 

E [^^J diAO) < \[\',\T\v.q)\\i^2 >=< pA;A'2|T|qAiA2 > . (6) 

Here 9 is the scattering angle between q and p and A = Ai — A2, A' = A'^^ — A2. The reduced 
rotation matrices dxy{6) are those defined by Rose [^. The matrix elements on the left 
hand side of Eq. (^) are independent of the third component of the isospin I3 and depend 
only on the absolute values p, q of the momenta. 

III. ON-SHELL SCATTERING 

The on-shell nucleon-nucleon cross section can be directly determined from the T-matrix 
amplitudes. In this case the extension to off-shell scattering is straightforward. Another 



possibility is to determine the on-shell phase shifts [21,30|. Doing so, an extension to the 



off-shell case is, however, unclear. Furthermore, the definition of the vacuum phase shifts, 
see e.g. refs. [0,^, has to be modified inside the medium to account for the modified 
unitarity relations. To be more precise, the definition of in-medium phase shifts should 
include the Pauli operator as pointed out in |[T6|jr7[l . To avoid such problems we will directly 
determine the cross sections from the matrix elements. The squared matrix elements are 
given as 



|T(p,g,^)r = EA 



2J + 1 
Air 



2J+ 1 



diy{9)ImTf{p,q) 



(7) 



The weighting factors /3j = 2 , i = l.A and jS^ = (3q = 4 arise from the sum over all helicity 
states. Using the orthogonality relation for the rotation matrices 
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dcosie)diA0)di:AO) 



2J+1 



-6 



jj' 



one obtains 



dn\f{p,q,e)\^ 



i=l J 



2 J 



47r 



[ReT^{p,q)f +(lmT^{p,q) 



(8) 



(9) 



In the on-shell case p = q the differential cross section follows from the matrix elements by 
the standard expression 



da 



■\T{q,q, 



'dn . 



(10) 



A. Free cross section 

The predictions of the Bonn potentials for free NN cross sections have in detail been 
discussed by Li and Machleidt [|2T|. Throughout this work we apply the Bonn A potential 
and the results of for the vacuum case are reproduced with high accuracy for both. 



differential and total cross sections. To demonstrate this feature, the results from ||2T| are 
included in Figs.|| and |^ where the total neutron-proton (T/ = 0.5(T/'^~° + T/'^~^)) and 
proton-proton (T/ = T/'^~^) cross section are shown. As found in the pp cross section 
is in particular at low energies significantly smaller than the np cross section. It should. 



however, be noted that in the present work as well as in ref. |^T|] the pp cross sections are 
not Coulomb corrected. 
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B. In-medium cross section 



As already pointed out in refs. |]18ip!9| a trivial medium modification of the cross sections 
arises from the in-medium masses entering into the kinematical term {M*Y / s* in eq. (p!OD 
which is due to the normalisation of the relativistic spinor basis and the incoming flux. This 
phase space factor reduces the in- medium cross section by the order of (M*/M)^ at small 
momenta. 

Besides the fact that one deals with dressed quasi-particles the essential feature of the 
Bethe-Goldstone or Bethe-Salpeter equation, respectively, is the occurrence of the Pauli 
operator. This means that the vacuum relations which connect the phase shifts with the 
real reaction matrix R are modified by the Pauli operator. In a schematic notation^ the 



in-medium reaction matrix R is connected to the T- matrix by |TD|,^| 



R-f = mRQT (11) 

which leads to a modified optical theorem 

|f|2 = = (7rg)-V^^I • (12) 

With Q = 1 the vacuum expressions are recovered. The modification of the optical theorem 
by the presence of the medium, in particular the appearance of the inverse Pauli operator 
which compensates at momenta below the Fermi surface for the vanishing ImT, has been 
discussed in []TB|,|T^,^ . It becomes clear from ([T^) that the use of the vacuum relations (with 



Q = 1) to extract phase shifts from the in- medium reaction matrix R is an approximation 
justified at low densities and/or high energies. In between the Pauli operator is essentially 
different from unity and appears in the denominator of eq. (|I2|). Neglecting here the 
influence of the Pauli operator will lead to an underestimation of the corresponding cross 
sections. The real and imaginary part of the T-matrix are related to the reaction matrix by 



For the clarity of the notation we suppress in eqs. (|llHl|) factors ^ which can be absorbed into 
R and T [^] and the (5-function originating from the principle value treatment of the Blankenbecler- 
Sugar, respectively the Thompson propagator. 
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We emphasise this point because, as will be seen in the following, we find the in-medium cross 
sections to be substantially less suppressed at low momenta than found by Li and Machleidt 



2T| whereas we obtain a good agreement with their results at high momenta. The reason 



for the deviations can be traced back to the different procedures used to determine the cross 
sections. As discussed in sec. 2 the squared matrix elements (|^) provide an unambiguous 
and direct method to extract the cross section. To determine phase shifts first, has to be 
done with caution since the occurrence of the Pauli operator modifies the corresponding 
phase shift relations in the medium ||T^. If neglected, as done in the approximation used in 



Ref. pT[, this effect leads in particular at low momenta to an underprediction of the cross 
section. To illustrate this effect in Fig. 4 we investigate the influence of the Pauli operator 
on the in-medium np cross section. For a fair comparison the density {kp = 1.4 fm~^) as well 
as the value of M* are chosen as in Ref. [^. One curve in Fig.|l] is obtained by switching 
off the Pauli operator in the Thompson equation (|^), i.e. setting Q = 1. It is clearly seen 
that the influence of the Pauli operator leads even to an enhancement of the cross section 
at momenta below ^180 MeV compared to the Q = 1 case which is due to the occurrence 
of Q in the denominator in Eqs. (12, 1^). As expected, at very small momenta the presence 
of the Pauli operator leads to a signiflcant suppression of the cross section. One should, 
however, keep in mind that the Pauli operator acts here only on the intermediate states in 
the Thompson equation ([l|) and not on the flnal states. Thus full Pauli blocking Q = 
reduces the full T-matrix to its Born part V. In the transport approach, on the other hand 
side, the Pauli blocking prevents also the scattering into occupied flnal states. Thus Fig. 4 
refers only to the Pauli effects in the intermediate states, but demonstrates the importance 
to account properly for the Pauli operator in the in-medium optical theorem (|12j). 
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FIG. 1. Total np in-medium cross section at kp = 1-4 fm ^. The result of the full calculation is 
compared to a calculation where the Pauli operator in the Thompson-equation for the intermediate 



scattering states has been switched off ((5 = 1). Also the corresponding result of Ref. ||2^ is shown. 



In the following we consider the in-medium cross sections at four different Fermi momenta 
kp = 1.1, 1.34, 1.7 and 1.9 MeV which corresponds to densities p = 0.090, 0.1625, 0.332 
and 0.4633 fm~^. For simplicity we denote these densities in the text as 0.5/1/2/3 po 
although these values do not exactly correspond to multiples of po = 0.1625 fm^^. Again 
all calculations are performed using the Bonn A potential. The in-medium mass M* entering 
into eq. (p!OD has in our calculation the values M* = 766.6, 646.7, 433.6, and 310.1 MeV. 
These values are slightly larger than those of ref. PTIITU]. Thus we expect also slightly larger 



values for the in-medium cross section, in particular at low momenta which are due to higher 
values for the kinematical factor (M*)^/s*. The reason for the different effective masses lies 
in different solution techniques for the Thompson equation (|I]). As discussed in Sec. 2, 
we apply a refined projection scheme in order to transform the T matrix from the two- 
particle cm. frame to the nuclear matter rest frame where the self-energy components 
are determined [|1^,|1^]. In the medium the on-shell T-matrix (|I|) depends on three variables, 
the relative momentum q of the initial states, the scattering angle 6 and the centre-of-mass 
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momentum P of the two-particle cm. frame relative to the nuclear matter rest frame. As in 
refs. PT],[TB|,0 we consider only the case where the two-particle cm. frame and the nuclear 
matter rest frame coincide, i.e. P = 0. 

In Fig. the differential np cross section at the different densities is shown at fixed 
relative momentum |q|=342 MeV which in the vacuum corresponds to a laboratory energy 
of i?iab = ^=250 MeV. The vacuum definition of -Eiab was used in [^] to compare the 
cross sections at different densities. The presence of the medium tends to make the np 
differential cross section more isotropic. At backward angles the cross section are decreasing 
with density. At forward angles the behaviour is more complicated: At moderate densities 
(p = 0.5/1 po) the cross section is reduced but at high densities (p = 2 and 3 po) a strong 
enhancement of the forward scattering amplitude can be observed. It is worth to notice 
that in this energy range we are in good agreement with the results obtained by Li and 
Machleidt [^. Similar results have been obtained at 0.5, 1 and 2po and also at 2po the 
cross section was found to be enhanced at forward angles compared to po. Going higher in 
density (3po) we find this effect even more pronounced. While the cross section stays now 
almost constant at backward angles it is strongly enhanced at forward angles. However, at 
3po the cross section is highly anisotropic and dominated by a p-wave component. Here we 
see a suppression of higher partial waves with increasing density. At 3po one needs partial 
waves up to at least J < 6 to approximate the full result (J = 12), at 2po the partial waves 
J < 4 are almost sufficient and at 3po the behaviour is like an s +p-wave with J < 1. Fig.|^ 
shows the same for pp scattering. Again our results are in good agreement with the findings 



of ref. In the 7 = 1 channel (pp) the cross sections are generally more isotropic than 
in the np channel. With rising density the cross sections are first decreasing (0.5po, Po) 
and then increasing. At 3po we observe a dramatic increase of the cross section at forward 
angles. In the pp cross section the contributions of higher partial waves are reduced with 
growing density, e.g. at po partial waves up to J < 4 reproduced the full result quite well 
whereas at 3po only contributions from J < 3 are relevant. 
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FIG. 2. Differential np in-medium cross section at fixed relative cm. momentum |q|=342 MeV 
= 250 MeV) at various densities. Experimental data from [^] for the free scattering are 
included. 
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FIG. 3. Differential pp in-medium cross section at fixed relative cm. momentum |q|=342 MeV 
= 250 MeV) at various densities. 

The peculiar behaviour at 3po seen in Figs. @, |^ can be understood from the presence of 
the mean field. In Figs. ^ § we investigated the density dependence of the cross sections at 
an equivalent relative cm. momentum q. This does, however, not correspond to equivalent 
energies. At finite density the laboratory energy Eiah{\(l\,p) = E{\q\,p) — M is given by 



and is therefore strongly modified by the presence of the mean field. At high densities the 
energy scale is stretched by the decreasing effective mass M*. This effect is responsible for 
the suppression of higher partial wave contribution to the differential cross section above 
2po if one compares the cross sections at identical cm. momenta but at essentially different 
incident energies. 

To illustrate this effect in fig.^ we show the density dependence of the differential np 
cross section at the same laboratory energy i^iab — 250 MeV. At comparable energies rather 
than comparable cm. momenta the difference in the differential cross section at po and 3po 
is now much less pronounced. 
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FIG. 4. Differential np in-medium cross section at laboratory energy i^iabCp) — 250 MeV at 
various densities. 



The suppression of the in-medium cross section at forward angles which occurs at higher 
densities can be understood from Fig.^ and eq. (p^. This figure illustrates the influence of 
the Pauli operator and the imaginary part of the T-matrix. It is seen that at po the imaginary 
part of T which contributes in the vacuum by about 50% to the forward scattering amplitude 
{6 = 0) is now strongly suppressed by the Pauli operator. This effects is maximal at low 
momenta and high densities. When q lies below the Fermi surface the imaginary part of T 
vanishes completely. Thus, Pauli blocking in the intermediate states makes the cross section 
more isotropic. 
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FIG. 5. Differential np in-medium cross section at fixed relative cm. momentum |q|=342 MeV 
in the vacuum (upper curves) and at kp = 1.34 (po) (lower curves). The contributions from the 
real part of the T-matrix are shown separately. 

In Figs. ^ and |^ we show the total in-medium np and pp cross sections in the considered 



density range 0.5po — 3po as a function of -Eiab, Eq. ([I^). Using this quantity the scale is 



considerably stretched compared to the vacuum expression ^ (used in [^). There are two 
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major aspects to be noticed: At high energies i^iab > 200 MeV we find good agreement with 
the previous calculations of ref . |^ . For np as well as pp scattering the cross sections reach 
asymptotic values around 15-20 mb. At high densities the cross section has the tendency 
to rise again with increasing laboratory energy. This behaviour is even more pronounced in 
the proton-proton channel and has been observed by other groups as well [ p!^ , pT| ] . 

At low energies the in-medium cross sections are considerably less suppressed than ob- 



served in ref. One reason are our somewhat larger values for the in- medium mass, but 
this effect is not sufficient to explain the deviations at low energies. As illustrated by Fig.|I| 
the differences can be understood by the low density/high momentum approximation made 
in which neglects the influence of the Pauli operator in the optical theorem. Besides 
this point, at 0.5po we see a small additional enhancement of the np cross section around an 
energy of 15 MeV which not present in the pp channel. A much stronger enhancement of the 
cross section at low densities has been observed by Aim et al. [l7|. In the flnite temperature 
approach of |]T^ this critical enhancement is attributed to the onset of superfluidity. Crucial 
for such a superfluid state are contributions from hole-hole scattering in the Pauli operator 
which are absent in the standard Brueckner approach (used here). However, as discussed 
in a signature of a bound pair state can appear at low densities even when hole-hole 
scattering is neglected in the Pauli operator. In the present calculations such an resonance 
like enhancement of the cross section is only seen in the / = amplitudes which correspond 
to the quantum numbers of the deuteron, i.e. the ^5*1, ^Di and the ^Si-^Di transition chan- 
nels. Therefore the low density enhancement of the np cross section can be interpreted as a 
precursor of a superfluid state and supports the flndings of [p5|] . 
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FIG. 6. Total np in-medium cross section at various densities as a function of i^'iab(p)- In addition 
the results of (squares) for the free cross section are shown. 
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FIG. 7. Total pp in-medium cross section at various densities as a function of -Eiablp)- I^i addition 
the results of (squares) for the free cross section are shown. 
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IV. HALF OFF-SHELL SCATTERING 



A. Matrix elements 

In the following we consider the case of half off-shell scattering. This means that the 
initial states with cm. momenta ±q are on their mass shell E*{q) = Vq^ + M*'^. However, 
now the momenta of the final states ±p can vary independently, i.e. |p| ^ |q|. Since the 
requirement of energy conservation fixes the final energies = {po, ±p} to po = ^^/s* = 
E*{q), these states are off energy shell. To obtain an impression how the off-shellness affects 
the T-matrix amplitudes we consider first the matrix elements. In the on-shell case the cross 
sections follow from the squared matrix elements \T{q, q) p by Eq. (10). The squared matrix 
elements |T(p, g)p = / \T{p,q,9)\'^ dQ are also those quantities which enter directly into a 
generalised transport equation as transition amplitudes for off-shell scattering |2^. In that 



case it makes more sense to speak in terms of transition amplitudes than in terms of cross 
sections. The latter ones are obtained from the transition amplitudes by integration over the 
final state spectral distributions. Thus, the cross sections depend crucially on the spectral 
width of the particles whereas the transition amplitudes themselves are independent of the 
spectral functions. 

In Fig.§ we show the iso-spin averaged matrix elements |T(p, g)p for A^A^ scattering as 
a function of p and q at nuclear matter densities 0.5/1/2/3 po- The amplitudes |Tp are 
given in fm^ and have to be multiplied by the factor (27r)^ (due to our normalisation of the 
T-matrix) when they are inserted into equation (0) to obtain cross sections. The diagonal 
{p = q) on-shell elements correspond to the total cross section. First of all, it is seen that 
the amplitudes |Tp show the same behaviour as the corresponding (on-shell) cross sections, 
namely a general decrease with momentum and density. In the cross section this tendency 
is just enhanced by the kinematical factor M*^/s* which decreases with density and also 
with momentum. 

The off-shellness of the final states now given by 

A^ = E*(q)-E*(p)^^l^ (15) 

which is maximal perpendicular to the diagonal. It can be seen from Fig.§ that at low 
momenta the matrix elements are only weakly affected by the fact that the outgoing states 



are off energy-sfiell. Tliis does not mean tliat tlie matrix elements do not cliange going away 
from tlie on-sliell point. Indeed, tlie variation of tlie matrix elements is considerable as will 
become even more clear further on. However, the dependence of the matrix elements is 
nearly symmetric in q and p, i.e. 

|T(p,g)|2^|T(g,p)|2 (16) 

and thus not strongly affected by the off-shellness of the outgoing states. In particular 
in the low momentum region the the matrix elements fall off symmetrically with increasing 
momenta p and/or q. Only at the highest density of 3po the asymmetry is rather pronounced. 
The off-shell variation of the outgoing states leads here to a resonance structure around 
q = 250 MeV where the amplitudes increase with different strength in p and q directions. 
The off-shell behaviour is similar in the np and pp channels, except that the resonance 
structure at 3po is more pronounced in the latter case. The somewhat stronger off-shell 
dependence of the pp scattering at high densities is also reflected in the stronger reduction 
of the corresponding total cross section shown in Figs. [TT|JT^. 



19 




FIG. 8. Squared iso-spin averaged half off-shell T-matrix elements |rp in [fm^] at various den- 
sities. The initial states with relative cm. momentum q are on-shell, the final states with relative 
cm. momentum p are off-shell. 

In order to obtain a more quantitative impression of the off-shell dependence next we 
investigate the deviation of the amplitudes 5|Tp from their on-shell values as a function of 
the off-shellness 5u} of the final states. 5|Tp is thereby defined as the relative deviation, i.e. 

|rfa,MP-|rto.fo = o)P 
^1^1 = |rfe,fo = o)|^ • 

For a better comparison of the different densities the variable 5uj, i.e. the energy shift of the 
final states with momenta ±p relative to their on-shell energies, is scaled by the effective 
mass 

<5a;=^(i?*(q)-E*(p)) . (18) 
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The range for the variation of Suj is constrained by the kinematical hmits of our calculations, 
i.e. < p, g < 400 MeV. The symmetric fall off with p and q in the low momentum range 
which is reflected in Fig.^ implies already the following tendency: For p < g, i.e. 5uj > 0, one 
expects an enhancement of the amplitudes, i.e. 5|Tp > 0, whereas for p > q, i.e. 6uj < 0, 
one expects a reduction of the amplitudes. For small values of q this behaviour is clearly 
seen from Fig.^. There the deviation of the isospin averaged matrix elements 5|Tp from the 
on-shell point is shown at fixed values of the on-shell momentum q = 100, 200, 300 and 400 
MeV. First of all we see that the variation of the amplitudes is pronounced. Within the range 
of 6uj = ±50 MeV the matrix elements can easily vary by more than ±100%. Secondly, we 
find that the pattern are similar at moderate densities 0.5po and po but become essentially 
different at large densities p > 2po- Thus, the off-shell behaviour reflects a strong medium 
dependence. The systematics of the different pattern is, however, quite complex. 

Let us first consider moderate nuclear matter densities. At low momenta q the amplitudes 
show an extremely strong variation around the on-shell point which reflects their steep and 
symmetric fall off (see Fig.§). As already mentioned, this results in a strong suppression 
at negative Su and an equally strong enhancement at positive 6uj. In the high momentum 
region the variation of the amplitudes is much weaker which results in a smoother and less 
pronounced dependence on 6uj. With increasing density the situation changes drastically 
and is even reversed at large q: now the amplitudes are strongly enhanced in the negative 6u} 
region and reduced at positive 600. This reflects the asymmetry of |Tp around the on-shell 
diagonal seen in Fig.H in the high momentum range. 

To summarise: The half-off- shell matrix elements show a pronounced density dependence. 
At moderate densities, however, the dependence on the incident on-shell momenta q and the 
final momenta p which are off-shell is to large extent symmetric in p and q. This implies that 
the matrix elements are not too much affected by the shift to off-shell energies but are mainly 
determined by the absolute values of the momentum states p and q. As a consequence, the 
off-shell matrix elements can be approximated with an accuracy of about 10 — 30% by the 
on shell points in the following way 



which follows from the symmetry assumption (piBI) . At large nuclear matter densities this 





(19) 
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symmetry is more strongly violated, i.e. by about 20 — 40% at 2po- Nevertheless, in view of 
the extremely large variation of the matrix elements with 5uj, ranging from almost complete 
suppression to an enhancement of more than a factor of two, the accuracy of this symmetry 
assumption is quite good in the considered 5uj interval. It can be applied in a straightforward 
way, e.g. in transport calculations and requires only the knowledge of the on-shell matrix 
elements. However, at 3po the amplitudes are highly asymmetric in p and q and thus 
approximation (|T9|) does no more hold. Here an accurate description requires the knowledge 
of the exact matrix elements. 

Fig.0 illustrates the validity of the symmetry assumption (|19D and the strength of the 
explicit dependence of the matrix elements on the energy shift 5uj. There the relative 
deviations 



.|^|2 _ |T(g,g)p-|T(p,g)p 
from the exact results |T(j9, g)p are shown as functions of 5uj, again for fixed values of q. 
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FIG. 9. Relative deviation of the isospin averaged T-matrix elements |Tp from their values at 
the on-shell point as a function of the off-shellness 5uj of the final states. 5|Tp is shown for various 
densities at fixed incoming relative momenta q = 100, 200, 300, 400 MeV. 
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FIG. 10. Relative deviation of the isospin averaged approximated T-matrix elements |Tp, Eq. 
([l^), from their exact off-shell values as a function of the off-shellness 6u! of the final states. (5|Tp 
is shown for various densities at fixed incoming relative momenta of g = 100, 200, 300, 400 MeV. 



B. Cross sections 

A cross section is obtained from the transition amplitudes |Tp by the division through 
the incoming flux 



Vl2 



F*{qi,q2) 



*4 



E*(qi)E*(q2) E*(qi)^*(q2) 
and the multiplication with the final state phase space factors 



(21) 
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da = -^(2vr)^5^(g, + q, - - P2mpm,qiq2)\'^,A{p,)^^A{p,) . (22) 

The 5^-function ensures energy-momentum conservation. In the general case where the final 
states pi,p2 are off-shell, A represents the full positive energy spectral function. In the 
quasi-particle approximation A reduces to the on-shell condition 

A{p) = 2n6{p*^ - M*^)2Q{p*) . (23) 

Thus the spectral function fulfils the sum rule 

^Mp) = ^ ■ (24) 
(27r) E*{p) ^ ' 

Here we choose spectral functions of a Breit-Wigner form 

A{p) = 2'R- '^-^ 2 • (25) 

Thus A satisfies the sum rule and in the zero width limit F i — > the quasiparticle approx- 
imation (|23| ) is recovered. In the two-particle cm. frame the integral (|22|) can easily be 
evaluated. The total cross section is then given by 
2M*^ 



z V r 
da= , --dn p^dp\T{p,q,e)\^ 

dp, (2E*(q)-po) Po 



((2E*(q) - po)2 - (p2 + M*2))^ + (2E*(q) - p^f [{pi - p2 _ M*2)2 + p2 ^2] 

which reduces to expression (|1^) in the zero width limit. As discussed in the width F is 
determined by the imaginary part of the self-energy 

M* 

r = ~ ^*(p) ^"^^^(^' p) + ^"^^o(p, p) (27) 

and depends on density and momentum. In the present approach F follows from the imag- 
inary part of the T-matrix. It can, however, only serve as an estimate for the full particle 
width inside the medium. As discussed in Sec. II the Bethe-Salpeter equation is solved in the 
quasi-particle approximation. For the determination of the T matrix this treatment is justi- 
fied since the condition F <^ E^ .^ = E* — J^q is readily fulfilled. However, since the standard 
Brueckner-Hartree-Fock approach accounts only for the particle-particle correlations of the 
Brueckner hole-line expansion JmS vanishes for momenta below the Fermi surface due to 
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(26) 



Pauli blocking. Long-range correlations which are usually treated in RPA-type approaches 
by a ladder summation in the particle-hole channel contribute to the spectral width 



but are not taken into account in the standard Brueckner-Hartree-Fock. RPA correlations 
lead further to a depletion just below the Fermi surface and an occupation of states above 
the Fermi surface. Thus, the DBHF results for F do not represent the full width but give a 
reliable estimate only at momenta well above the Fermi surface where particle-particle cor- 
relations can be regarded as the the dominant contributions to F. In the considered density 
and energy range F ranges from about 10 MeV at 0.5po to more than 40 MeV at 3po as the 
outcome of the present DBHF calculations [0. 

Thus, the cross sections shown in Figs. |Tl| (np) and |1^ (pp) give an estimate for the 
off-shell dependence of the total cross section. For this purpose we choose two typical values 
for F which cover the range of the spectral width in nuclear matter as it is predicted by 
DBHF calculations, i.e. F = 10 MeV and 40 MeV. The present results are obtained with a 
constant, i.e. momentum independent F. For a better comparison of the various densities 
the cross sections are shown as a function of 2q^/M which corresponds in the vacuum to 
the laboratory energy. It can be seen from Figs. [11], |12] that the off-shell dependence of the 
total cross section is moderate. The averaging over the Breit-Wigner distributions leads to 
significantly smaller off-shell effects than are seen in the scattering amplitudes. Compared 
to the on-shell value the cross sections are generally reduced, in the case of F = 40 MeV by 
about 15-20%. At small nuclear densities this reduction is most pronounced at low momenta 
whereas but at higher momenta there is no sizable effect. At large densities the reduction 
is more pronounced at high momenta. Only for large values of F ~ 100 MeV (not shown 
here) a reduction of the cross section of more than 50% can be reached. Differences in 
the off-shell structure of the NN interaction, e.g. using BONN C, do not lead to essentially 
different results. This is consistent with the observations made in Ref. |2T| where it has been 
demonstrated that the on-shell cross sections are not much affected by the use of different 
parameterisations like BONN B or C. 
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FIG. 11. Total np cross section at various densities function of 2q^/M. The quasi-particle 
approximation (F = 0) is compared to the case where the final states have a finite spectral width 
of F = 10, 40 MeV. 
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FIG. 12. Total pp cross section at various densities as a function of 2q^/M. The quasi-particle 
approximation (T = 0) is compared to the case where the final states have a finite spectral with of 
r = 10, 40 MeV. 



V. SUMMARY 



In the present work we investigated the in-medium nucleon-nucleon cross section within 
the relativistic (Dirac) Brueckner-Hartree-Fock approach. We considered both, on-shell 
scattering and the more general case where the final momenta are allowed to be off energy 
shell. The in-medium cross sections can serve as input for transport calculation of heavy ion 
collisions. Information on the off-shell dependence of the scattering amplitudes is required 
when one intends to go beyond the quasi-particle approximation in order to obtain a more 
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realistic description of transport phenomena for particles with finite width. In a dense 
hadronic environment created, e.g., in heavy ion collisions also 'stable particles' like nucleons 
acquire a finite spectral width. Throughout this work we applied the BONN A potential 
as the nucleon-nucleon interaction since this potential yields the most reasonable saturation 
properties for nuclear matter. 

Concerning the on-shell scattering we find a qualitative agreement with previous inves- 



tigations of Li and Machleidt however, less reduced in-medium cross sections at low 
energies. The reason therefore lies in the fact that we account for modifications of the optical 
theorem due to the presence of the medium which have been neglected in the approximation 
made in In the np cross section an additional low density enhancement appears which 
can be interpreted as the precursor of a superfluid state. The present approach was then 
extended to the case where incoming and outgoing momenta of the scattered nucleons can 
vary independently and - due to energy-momentum conservation - the final states are off en- 
ergy shell. The resulting T-matrix elements or transition amplitudes show a strong variation 
around the on-shell point. The shape of the transition amplitudes depends, however, mainly 
in a symmetric way on the incoming/outgoing momenta. The fact that the final states are 
off-shell, the incoming ones on-shell, plays thereby a minor role. This allows to approximate 
the half-off-shell matrix elements in a suitable way by their on-shell values. Such an approx- 
imation works well at moderate nuclear matter densities but at large densities the precise 
knowledge of the amplitudes is required. In a conservative estimate finite width effects lead 
to a reduction of the total cross sections by about 20% compared to the on-shell scattering. 
Off-shell effects are more pronounced at larger nuclear densities. 
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